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Abstract 

We present the three loop contribution (order e 4 ) to the pressure of massless quan- 
tum electrodynamics at nonzero temperature. The calculation is performed within the 
imaginary time formalism. Dimensional regularization is used to handle the usual, 
intermediate stage, ultraviolet and infrared singularities, and also to prevent over- 
counting of diagrams during resummation. 
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The equation of state (EOS) of relativistic quantum electrodynamics (QED) at 
nonzero temperature (T) and chemical potential (fi) is of relevance in several astro- 
physical contexts |], 0, |, It was obtained by Akhiezer and Peletminskii to the 
third order (e 3 ) more than three decades ago while the fourth order (e 4 ) contribution 
at T = 0, but nonzero was then obtained for QED and quantum chromodynam- 
ics (QCD) by Freedman and McLerran 0, and Baluni ||. However, to date, the 
analogous 3-loop calculation for T ^ has, to our knowledge, not been performed, 
presumably because of the greater technical difficulty in dealing with overlapping 
diagrams in the presence of Bose-Einstein and Fermi-Dirac statistical factors. 

In this Letter we present the order e 4 contribution to the pressure of a QED plasma 
at T ^ (and fi = 0), but for the case of massless electrons. This will also be the 
leading contribution at high temperature for realistic QED with massive electrons. 

Our motivation for this high order calculation is two-fold. Firstly, for phenomeno- 
logical applications it is helpful to know how big corrections to the lower order EOS 
can be. Secondly, it serves as a prototype to illustrate techniques which may be be 
used to perform similar calculations in QCD. Recall that in QCD asymptotic free- 
dom suggests [H that at high T and/or density hadronic matter will transform into 
a weakly interacting quark-gluon plasma, a novel state of matter that is currently 
under intense study ||. The EOS for the high-temperature phase of QCD was de- 
termined by Kapusta to third order (g 3 ) [[jj while Toimela || has extracted also the 
gHn g piece. However, the normalization of the logarithm in the last term requires 
knowledge of the 3-loop contribution which is still lacking. 

Returning to QED, let us introduce our notation and conventions. We employ the 
imaginary time formalism whereby the energies take on discrete Matsubara values, 
go = innT, n being an even (odd) integer for bosons (fermions). Dimensional regular- 
isation is be used to handle the ultraviolet (UV) and infrared (IR) singularities which 
occur at intermediate steps. The D dimensional vector, = (q , q) is contracted 
with a Minkowski metric, Q 2 = q$ — (f. In order to keep track of the even (odd) 
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Matsubara frequencies, we introduce the following notation: 



f[dq]=T £ / 



d D ~ l q 



qo,even ' 



(2 



7T 



(27 

qo,oaa \ 

The fermions are kept as four- component objects, Tr(7 M 7„) = 4g^, and for simplicity 
we will work in the Feynman gauge so that the gauge propagator is g^/K 2 . Renor- 
malization via minimal subtraction ensures that the coupling constant is gauge-fixing 
independent, and hence so will then be our final answer for the pressure |3[]. 

Before discussing the 3-loop calculation, let us summarize the lower order results 
for the pressure of QED with N massless Dirac fermions : 

P = P + P 2 + P 3 + 0(e 4 ), (1) 



where 



^(l + IiV), (2) 



P = 5e2T4iV (3) 
2 288 ' 1 ; 
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The ideal gas contribution P is determined by the one- loop diagrams in Fig. 1. 
The first correction, P 2 , is given by Fig. 2. The nonanalytic (e 2 ) 3//2 contribution, 
P3, is a consequence of Debye screening [^]. In a perturbative expansion using bare 
propagators one discovers power-like IR singularities in diagrams such as Fig. 3, cor- 
responding to the n = Matsubara frequency of the photon propagator. The di- 
agrams in Fig. 3 are singular because the electric polarization operator behaves as 



Hoo^O) Q — > 0) = m 2 , where m is the electric screening mass to lowest order. Summing 
the infrared divergent pieces of the diagrams in Fig. 3 yields 

T f d 3 q f, m 2 N m 2 ^ 



2 J ( 2.)3V ln(1 + ^ ) -^J' (6) 
Although (H) is UV and IR finite and may be evaluated directly to give @, it 
will be helpful to reconsider it using dimensional continuation. Then |10| the second 
term in (|6]) vanishes and the first term gives 

T r (l-D\ « ( " 1V2 



-r, 
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which is P3 as D — > 4. Note that in dimensional regularization (DR) each single 
integral of the series in ([5]) vanishes. Neverthless the result in (§) is physically correct, 
and mathematically nonzero, because the sum in (|J) must be done before the integral. 
This is obvious once one starts from an expansion of the pressure in terms in terms 
of the full propagator 0, |, |J. Though the infinite sum must be performed first, we 
may separate out any finite number of terms from the sum, and for these we may 
use DR to deduce a zero contribution. This explains why the second term in ([]) may 
be dropped in DR and we will exploit this fact further below. 

Consider now the 3-loop diagrams. The order e 4 N diagrams are shown in Fig. 4 
and 5. However, for massless fermions, the Ward identity Z\ = Z2 implies the mutual 
cancellation of the counterterm diagrams. Thus the sum G\ + G2 (Fig. 4) is UV 
finite. After performing the spinor traces, some algebraic manipulation, and the use 
of scaling arguments such as in [TT|, we obtain 

? = ^(2(l-2<"-->) * I + (20-3D)*) 

+ (D - 2) ! ( 2 i/ 3 - / 2 (/, - M 2 ) . 

Here \i is the mass parameter of DR and e is the dimensionless, renormalized coupling. 
We have scaled all the momenta by 1/T so that the integrals are dimensionless (i.e. 



T = 1 there) and are defined by 
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The integral Hi occurs in the 3- loop evaluation of the pressure in c/> 4 theory ||12|| , 
and H 2 , being simply the fermionic analog of Hi, may be analysed in a similar 
manner. The only new integral left is if 3. Let us however first discuss the order 
e 4 N 2 diagrams. The UV singularity of G3 (Fig. 6a) is cancelled by the photon wave- 
function renormalization through diagram, X\ (Fig. 6b). Diagram G3 also has an 
IR singularity which is precisely the first term of the series in eq.(^). Since this 
term has already been considered there, it should be subtracted from G 3 to avoid 
overcounting. However, as discussed earlier, this single term by itself vanishes in DR, 
so double-counting is automatically avoided. We have 

G 3 = ^ ji3D-8 ^-2D 16 ( {D _ 4) ^ f 2 + P_4) ^ +4H \ (g) 



and 

Xi = -(Z 3 -l)e 2 N{D- 2) T 3D ~ 8 ^~ 2D h (2bi - h) 



'T\ 4 ~ D 



where 

e 2 N 

3 ~ 1 ~ 6tt 2 (D-4) 

to leading order. The factor ( — J will generate a e In Tj \i term on expansion, but 
this will be reabsorbed later in the definition of the temperature dependent coupling 
e(T). The new integral in @ is 
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We now sketch our evaluation of H 3 fl7|). The frequency sums are first rewritten 
in terms of contour-integrals as in |2|, |3| but we do not separate the T = and T^O 
parts before going to the "phase-space" representation. In this way we obtain 



# 3 = J i + K x + L 



i- 



The piece L\ contains integrals which can be performed analytically while K\ is an 
integral similar to H\. The difficulty lies in 



Ji=f ^|^(|^ ^(^ 2 )^(^ 2 )^(Q 2 )^ n go (iV Po +n P0 ) J2 (-^,7), 

V / CT,7=±1 

with 

{CT,1> K ■ Q K ■ Q + P ■ (aK + 7Q) ' 

and 

The integral J\ is UV finite but has a collinear singularity. We use the method of 
Sudakov decomposition, a technique which is well known at T = (see for example 
||13|| ) but appears to be novel in this context, to extract the pole and the finite part 



of this integral. The details are lengthy and will be presented elsewhere [JH|. Here 
we only state the final result, 

Jl = 128^ ((£>-4) + Cl ) ' 

where 

n = -0.7167667897 ± 10" 10 , 
ci = 3.936 ±5 x 1(T 3 . 

The residue r± is obtained as a finite two-dimensional integral which we have then 
evaluated numerically to high precision (relative error 10™ 10 ), and the pole in J\ above 
cancels with the other poles contributing to the sum of G\ + G<i (again to the same 
precision). For if 4 , the analysis is similar to H 3 but more tedious because of the 
doubled propagator 1/(Q 2 ) 2 . 

Collecting all the pieces, we obtain 

P 4 (gi + G 2 + G 3 + XQ 

e 4 T 4 e 4 T 4 

= ^(0.4056) - *M«L + , 5 hlV (9) 
vr 6 v 7 \ 7T 6 6tt 2 x 288 /i / v ; 

The pressure up to and including order e 4 then follows from eqs. (|I]-fi|) and @. It 

may be rewritten in terms of the (one-loop) renormalization group invariant coupling, 

at the energy scale T, given by 

Defining a (T) = e 2 (T) / 4ir we finally arrive at 

P_ ^ 7 Svr 2 q(T)iV 2vr 2 / a(T)iV \ 3/2 

T 4 45 1 4 ' 72 vr 9^ V ^ / 



+ ( 0.658±0.006 _ 757±0Q04 ^ /™ 2 + 0(a(T)5/2 



(10) 



This (|H]) is our expression for the three-loop pressure of QED with N electron 
flavours at high temperature; if m e is the (zero temperature) electron mass, we require 
m e /T <^ Oi{T) 2 so that mass corrections are subleading to the terms displayed in (|T0|). 
Real world QED corresponds to iV = 1 and in the regime where a{T) <C 1, the three 
loop contribution is found to be indeed a small correction. However, one should note 
that since perturbative QED is not asymptotically free, the effective coupling a{T) 
increases with temperature (albeit slowly) so that at sufficiently high temperatures the 
fourth order contribution becomes relevant. Defering further discussion and potential 
applications of (|T0D to a later stage |14j], we mention that the unknown e 5 contribution 
in eq. (|It]) is a higher order analog of the e 3 plasmon term and is also calculable |TJ) 



We conclude by summarizing the steps leading from the Feynman diagrams to 
the result (|T0D : (i) algebraic reduction of the integrals, (ii) evaluation of frequency 
sums by a contour-integral algorithm, (iii) evaluation of the final phase-space-like 
integrals, in particular using the method of Sudakov variables for UV finite integrals 
with a collinear singularity, and, (iv) use of dimensional regularization to simplify 
the prevention of double counting during resummation. We hope that the method- 
ology adopted here opens the way for a similar, and long awaited, calculation in QCD. 
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Figure Captions 

Fig-l : 

Contribution to the ideal gas pressure. The wavy line represents the photon propa- 
gator. 
Fig-2 : 

The two loop diagram. 
Fig-3 : 

Diagrams contributing to the e 3 plasmon term. The self-energy insertions are n 00 (0, 0) 
and the photon is static, q = 0. 
Fig.4 : 

The order e A N contributions : G 1 and G 2 . 
Fig-5 : 

Ultraviolet counterterm diagrams for Fig.4. 
Fig.6 : 

Fig. 6a is the e A N 2 contribution (G 3 ) while Fig.6b is the corresponding counterterm 
diagram X 1 . 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9405343v2 
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